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ABSTRACT 


This  thesis  deals  with  analytical  predictions  and  experimental  determination  of  the 
modal  frequencies  and  shapes  associated  with  the  arm  portion  of  the  Naval  Postgraduate 
School's  Flexible  Spacecraft  Simulator  (FSS).  A  description  of  how  piezoceramic  sensors 
and  actuators  are  incorporated  in  finite  element  modeling  is  presented.  A  MATLAB™ 
code  conducting  the  finite  element  modeling  of  the  arm  is  used  to  verify  the  modal 
frequencies  generated  via  Structural  Dynamic  Research  Corporation's  I-DEAS™ 
software.  Modal  testing  is  conducted  with  the  I-DEAS  Test  package  to  determine  the 
first  four  physical  modal  frequencies  and  shapes  so  that  we  may  compare  them  with  the 
analytical  results.  The  results  of  the  testing  indicated  that  finite  element  analysis  predicted 
modes  one  and  three  within  an  average  of  23.1%,  and  modes  two  and  four  within  an 
average  of  4.4%. 
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L  INTRODUCTION 


The  often  times  prohibitive  cost  for  putting  payloads  into  orbit  has  lead  the 
scientific  and  engineering  community  to  design  hghtweight  yet  dynamically  stable 
structures.  Accurately  predicting  the  natural  frequencies  of  a  structure,  corresponding 
mode  shapes  and  modal  damping  is  integral  to  the  design  and  dynamic  control  of  such 
stmctures.  Lightweight,  flexible  structures  inherently  possess  low  resonant  frequencies 
which  can  easily  be  stimulated  by  a  variety  of  spacebome  excitation  sources  such  as 
attitude  control  thrusters  and  reaction  wheels. 

Previous  work  done  on  flexible  structural  control  utilizing  Positive  Position 
Feedback  (PPF)  and  Strain  Rate  Feedback  (SRF)  by  Newman  and  Feuerstein,  [Refs.  1 
and  2]  mandated  an  accurate  knowledge  of  modal  characteristics  for  successful  control 
law  design.  This  previous  work  in  modal  control  provides  the  stimulation  for  this 
research. 

The  most  common  method  of  modal  frequency  and  shape  prediction  is  some  form 
of  the  Finite  Element  Analysis  (FEA)  approach.  This  method  exploits  digital  computing 
power  and  provides  a  formulation  for  the  analysis  of  dynamic  characteristics  of  flexible 
structures.  By  regarding  a  complex  structure  as  a  composite  of  finite  elements,  the 
method  requires  that  displacements  be  compatible  and  internal  forces  in  equihbrium  at 
points  known  as  nodes,  which  are  shared  by  several  elements.  The  result  being  the  total 
structure  should  act  as  if  one  continuous  entity. 

The  Naval  Postgraduate  School  Flexible  Spacecraft  Simulator's  arm  provides  a 
testbed  for  modal  analysis.  It  provides  a  reasonably  accurate  representation  of  a  very 
flexible  arm  such  as  that  which  might  be  utilized  for  a  deployed  antenna,  possibly  imposing 
very  accurate  pointing  constraints.  This  arm  employs  nearly  collocated  piezoceramic 
actuators  and  sensors  for  the  damping  and  control  of  vibration.  The  FEA  method  is  well 
documented,  [Refs.  3  and  4].  What  is  initially  explored  is  how  the  physical  properties  and 
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forces  inq)lied  to  a  structure  via  the  incorporation  of  these  piezoceranric  sensors  and 
actuators  modifies  the  FEA. 

Many  commercial  software  FEA  packages  are  available  to  assist  with  the 
predictions  described  above.  For  the  purpose  of  this  thesis,  the  modal  prediction  package 
provided  by  Structural  Dynamic  Research  Corporation's  I-DEAS  software  is  used.  A 
comparison  of  this  output  to  a  simple  MATLAB  code  which  utilizes  standard  FEA 
techniques  [Ref.  4  ]  is  then  undertaken. 

Finally,  modal  analysis  of  the  arm  is  conducted  experimentally  so  that  we  may 
determine  the  accmacy  of  analytical  predictions.  This  modal  testing  utilizes  the  testing 
package  incorporated  within  the  I-DEAS  software.  The  arm  is  excited  at  several  nodes 
with  an  impact  hammer  with  a  built-in  force  transducer.  Accelerations  of  the  arm  are 
measured  via  an  accelerometer  whose  signal  has  been  conditioned  and  amplified.  These 
analog  time  domain  signals  are  discretized  and  then  changed  into  the  fi'equency  domain  via 
the  discrete  Fourier  transform  The  analyzed  signals  can  then  be  manipulated  in  a  variety 
of  ways  to  produce  such  information  as  natural  fi'equencies,  mode  shapes  and  damping 
ratios. 
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n.  UNITE  ELEMENT  ANALYSIS  WITH  PIEZOCERAMIC 
ACTUATORS  AND  SENSORS 


A.  ANALYTICAL  BACKGROUND  AND  ASSUMPTIONS 

For  the  purpose  of  discussing  the  procedure  with  which  FEA  methods  are 

modified  to  take  into  effect  the  physical  and  active  characteristics  of  an  element 
incorporating  piezoceramics,  consider  the  co-located  sensor/actuator  pair  on  the  n-th 
element  of  an  isotropic  flexible  beam  as  depicted  in  Figure  1  below.  All  physical  and 
electrical  properties  and  dimensions  of  the  actuator  and  sensor  are  assumed  identical. 

The  neutral  axis  is  aUgned  with  the  centerline  of  the  beam  The  beam  is  Jfree  to  translate 
in  the  +/-  number  3  axis  direction  and  rotate  about  the  number  2  axis  only. 


Figure  1.  Piezoceramic  Actuator  and  Sensor,  n-th  Element 
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B.  GOVERNING  EQUATIONS 


1.  Passive  Contributions 

We  first  consider  a  structural  element  as  shown  in  Figure  2.  The  equations  of 

motion  for  the  element  are  of  the  form: 

Mq  +  Kq  =  ^  (1) 

where  M  and  A!"  are  the  elemental  mass  and  stiflfiiess  matrices,  respectively,  q  represents 
an  element  displacement  vector  of  the  form: 

W«-i 

0M-1 


where  w„-i  and  are  the  transverse  deflections  of  the  left  and  right  end  of  element  n, 
respectively,  and  6„-i  and  0„  are  the  rotations  of  the  left  and  right  end  of  element  n, 
respectively.  These  displacements  are  shown  in  Figure  2  below. 


Figure  2.  n-th  Element  Displacements 


By  deriving  the  equations  of  motion  for  the  element  using  Lagrange's  equations,  which 
reduces,  iu  turn,  to  determining  the  kinetic  energy,  the  potential  energy,  and  the  virtual 
work  expressions  in  terms  of  the  nodal  cordinates,  consistent  element  mass  and  stiflfiiess 
matrices  can  be  generated. 

The  consistent  element  mass  and  stiflfiiess  matrices  are  given  by  the  following: 
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156  22h  54  -\?>h  ' 

22h  4h^  13/2  -3/22 

"420  54  13/j  156  _22/2 

-13/2  -3/2^  -22/2  4/2^ 


12  6/2  -12  6/2  ' 

6/2  4/22  _6;j  2/22 

"  -12  -6/2  12  -6/2 

6/2  2/22  _6;j  4/j2  ^ 


where  m  is  the  hnear  mass  density  of  the  element,  h  is  element  length,  E  is  Young's 
modulus  for  the  element  and  I  is  the  element  moment  of  inertia  about  the  bending  axis. 
[Ref  4  ] 

In  order  to  accoimt  for  piezo  contribution  to  the  mass  matrix,  m  must  be  equal  to 
the  linear  mass  density  of  the  beam  plus  the  linear  mass  density  of  the  piezoceramic 


actuator  and  sensor,  or. 


To  account  for  piezo  contribution  to  the  stifBiess  matrix,  the  El  term  must  be  of  the  form; 

El  —  E beamJ-heam  "I"  E piezoJ piezo  (6) 

Since  the  neutral  axis  is  aligned  with  the  beam  centerline,  Iheam  is  simply: 

heam  =  U  (about  bending  axis)  (7) 

where  w  is  the  width  of  the  element  and  tb  is  the  thickness  of  the  structure. 

By  employing  the  parallel  axis  theorem  to  determine  1  piezo ,  we  arrive  at  the 
following: 

W  =  2[^  +  («"p)(t  +  t)'] 

where  tp  is  the  thickness  of  the  piezoceramic. 

After  simplification. 


Iniezo  —  2wt. 


±  ^ 

A  n 


Ebeam  and  E piezo  are  sin4)ly  the  Young's  modulus  values  for  the  beam  and  piezo, 
respectively. 
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Once  the  particular  elements  incorporating  the  piezos  are  modified  as  above,  all 
other  aspects  of  constructing  the  global  mass  and  stifiBiess  matrix  for  an  entire  structure 
incorporating  multiple  elements  are  straightforward  [Ref  4],  The  global  version  of 
Equation  (1)  will  now  present  an  eigenvalue  problem  whose  solution  will  reveal  modal 
frequencies  and  shapes. 

2.  Active  Contributions 
a.  Sensor 

Piezoceramic  actuator  and  sensor  theory  is  well  documented  [Ref  5], 
Voltage  output  from  the  piezoceramic  sensor  is  given  by  the  following  equation; 

Vs  =  tp(Epd3i/D)(ei  +  S2)  (10) 

where 

tp  is  the  thickness  of  the  piezoceramic  sensor 
Ep  is  Young's  modulus  for  the  sensor  in  N/nt  ^ 
dsi  is  the  lateral  strain  coefficient  in  nW  or  Coul/N 
D  is  Abs  Permittivity  in  Farad/m  or  N/V^ 

Si  and  82  are  longitudinal  and  lateral  strain,  respectively 

Since  longitudinal  strain  S 1  and  the  lateral  strain  82  are  related  by  Poison's  ratio,  V,  given 
by: 

SiV  =  -S2  (11) 

Equation  (10)  may  be  written  as: 

Vs  =  tp(Epd3ilD)il-v)si  (12) 

Figure  3  illustrates  the  sensor  mode  of  the  piezoceramic.  The  illustration  on  the  left  of 
Figure  3  depicts  the  case  when  the  sensor  is  subject  to  lateral  tension  and  vertical 
compression.  The  right  illustration  shows  the  sensor  in  vertical  tension  and  lateral 
compression.  Produced  voltage  polarities  under  these  conditions  for  the  indicated  poling 
direction  are  as  shown. 
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Figure  3.  Hezoceramic  Sensor 


In  generalized  coordinates,  sensor  voltage  produced  at  the  n-th  element  is 
represented  by  the  following,  [Ref  5]: 

f;  =  F,[  0 1 0 1  ]?  (13) 

where  Vs  is  given  by  Equation  (12)  and  the  displacement  vector  g,  is  given  by  equation 

(2). 

b.  Actuator 

The  equation  with  which  the  general  behavior  of  the  element,  exposed  to 
external  excitation  forces  and  also  the  forces  uiq)osed  by  the  actuator,  is  represented  by 
the  following: 

Mg  +  K.g  —  F external  "t"  F actuator  (14) 
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where  F actuator  is  given  by  the  following  expression,  [Ref  5]: 

F actuator  =  ~BV a  (15) 

Va  is  the  voltage  applied  to  the  actuator  and  5  is  a  4  x  1  vector  possessing  moment 
correction  terms  whose  units  are  in  Newton-  meters/volt: 

0 
bi 
0 
b4 

with 

bi  =  -d-i\Epw(-^  +  y)  and  b^  =  -bi  (17) 

where  aH  variables  are  as  in  Equation  (10)  and  w  is  the  actuator  width.  Figure  4 
illustrates  the  actuator  mode  of  the  piezoceramic.  The  left  side  illustration  demonstrates 
vertical  contraction  and  the  corrective  moment  producing  lateral  expansion  when  apphed 
voltage  is  of  opposite  polarity  fi'om  the  poling  voltage.  The  right  side  demonstrates 
vertical  expansion  and  the  corrective  moment  producing  lateral  contraction  when  a 
voltage  of  the  same  polarity  as  the  poling  voltage  is  apphed  to  the  piezo. 
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Figure  4.  Piezoceramic  Actuator 


c.  Coupling  Sensor  Output  to  Actuator  Input 
The  last  procedure  is  to  determine  the  control  scheme  which  wiU  couple 
the  measured  sensor  voltage  output  to  that  apphed  to  the  actuator  in  order  to  produce  a 
restoring  moment  to  the  element  and  ultimately  the  entire  structure  itself.  We  will 
examine  a  simple  implementation  of  velocity  feedback,  where 

V,  =  kV'.  (18) 

A:  is  a  voltage  gain  factor.  Combining  all  variables  in  Equations  (2),  (13)  and  (15)  will 

ultimately  lead  to  the  following  overall  e?q)ression  for  F actuator  • 

r  A  r  , 

0  M>„-l 

Factor  =  -k\  FJ  0  1  0  1  (19) 

b4  Q„ 
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With  the  variables  bi 
form: 


Z>4  and  Vs  as  before.  Equation  (19)  can  now  be  written  in  the 


actuator 


=  -Cq 


(20) 


where  C  is  a  4  x  4  matrix  given  by: 

C  =  -ki 


0 

bi 

0 

Z)4 


F.[  0  1  0  1  ' 


(21) 


With  no  external  excitation  forces  present,  Equation  (14)  may  now  be  rewritten  in 
the  general  damped  form: 

Mq  +  Cq  +  Kq  =  0  (22) 


In  the  construction  of  a  global  system  of  i  number  of  elements,  the  global  C 
matrix  will  be  null  in  all  segments  except  those  which  possess  actuators.  Therefore,  the 

passive  ri  ampin g  of  the  structure  is  neglected.  Equation  (22)  then  presents  a 
conventional  eigenvalue  problem  which  will  reveal  modal  frequencies,  modal  vectors  and 
darcping  ratios.  [Ref  6] 
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m.  FINITE  ELEMENT  ANALYSIS  COMPUTER  PREDICTIONS 


A.  I-DEAS™  FEA  MODELING 

I-DEAS  is  a  menu  driven,  Windows™  type,  software  package  developed  by 
Structural  Data  Research  Corporation.  Its  finite  element  analysis  and  testing  packages  are 
just  two  of  several  structural  design  tools  available  to  designers  and  engineers. 

The  first  step  in  conducting  FEA  with  I-DEAS  is  to  model  the  structure  under 
consideration.  Figure  5  depicts  the  mesh  geometry  for  the  Flexible  Spacecrafl;  Simulator's 
arm  used  for  the  analysis.  This  flexible  arm  rests  on  a  granite  table  which  simulates  a 
fiiction  fi'ee  environment  and  is  supported  at  the  elbow  and  tip  by  two  airpads  (nodes  10 
and  19).  Node  1  is  rigidly  fixed  to  the  table.  Mesh  geometry  is  such  that  the  arm  is  fi'ee 
to  translate  in  the  +/-  number  1  and  3  axis  and  rotate  about  the  munber  2  axis  only. 

The  arm  is  constructed  of  7075  T-6  aluminum,  with  a  cross  section  of  1/16"  x  1" 
(th  X  w),  and  has  several  mass  intensifiers,  which  are  treated  as  point  masses  for  the 
purpose  of  the  analysis,  spaced  throughout  the  beam  Elements  1,  10,  11,  and  18  have 

two  wafers  per  side  of  Navy  Type  11  PZT  bonded  to  the  beam  surface.  Table  1  below 
cites  the  material  and  electrical  properties  of  this  piezoceramic. 


Quantity 

Description 

Units 

Value 

d^i 

Lateral  strain  coefiF 

m/F  or  CoulIN 

1.8e-10 

Epiezo 

Young's  modulus 

Pascal 

6.3  elO 

V 

Poison's  ratio 

n/a 

0.35 

D 

Abs  Permittivity 

Farad! m  or  NiV^ 

1.5  e-8 

tp 

piezo  thickness 

m 

1.905  e-4 

Ppiezo 

Mass  density 

kgirrF 

7.7  e  3 

Table  1.  Navy  Type  11  PZT  Properties 
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ARM  STRUCTURE  ON  GRANITE  TABLE 
(Airpads  at  nodes  10  and  19) 


n  =  nodes 
e  =  elements 


n-2  n-3  n-4  n-5  n-6  n-7  n-8  n-9 


BEAM  MATE:  7075  T-6 

-  CROSS  SECTION:  1/16  "  x  1",  (tb  x  w) 

-  DENSITY:  0.003144  slugs/in3 
-E  =  1.029E7psi 


e-10 

n-11 

e-11 

n-12 


MASS  INTENSEFIERS: 


AIRPAD: 


=  0.0317  slugs 
=  0.0252  slugs 


AIRPAD: 


=  0.0567  slugs 


NODAL  COORD  (1,3  AXIS)  in. 


n-1 

(0,0) 

n-11 

(26.5,-2.5) 

n-2 

(4.75,0) 

n-12 

(26.5,-5.0) 

n-3 

(7.25,0) 

n-13 

(26.5,-9.87) 

n-4 

(9.75,0) 

n-14 

(26.5,-12.37) 

n-5 

(12.25,0) 

n-15 

(26.5,-14.87) 

n-6 

(14.75,0) 

n-16 

(26.5,-17.37) 

n-7 

(17.25,0) 

n-17 

(26.5,-19.87) 

n-8 

(19.75,0) 

n-18 

(26.5,-21.5) 

n-9 

(24.0,0) 

n-1 9 

(26.5,-26.5) 

n-10 

(26.5,0) 

•  ACCELEROMETER  =0.005  slugs 

NAVY  TYPE  II PZT  ON 
ELEMENTS  e-1,  e-10,  e-1 1,  e-18 


Figure  5.  Finite  Element  Analysis  Mesh  for 
Flexible  Spacecraft  Simulator  Arm 


J  n-13 
e-13 
n-14 

e-14 
I  n-15 

e-1 5 

3  n-16 

e-16 

D  n-17 
e-1 7 
1  n-18 


The  mesh  geometry  (all  nodal  coordinates),  beam  element  geometry,  element 
material  properties,  point  mass  locations  and  values,  and  structural  restraints  are  all  input 
into  the  I-DEAS  modeler.  Once  this  model  is  stored,  modal  analysis  can  be  conducted. 
I-DEAS  constructs  element  mass  and  stiflBaess  matrices,  assembles  the  global  mass  and 
stiEBiess  matrices  for  the  structure  and  then  calculates  the  modal  frequencies  and 
associated  mode  shapes.  Modal  analysis  was  conducted  on  models  with  and  without  the 
PZT  wafers  on  elements  1,  10,  11,  and  18.  Appendix  A  contains  the  I-DEAS  printouts 
for  the  model,  and  the  deformed  geometry  for  the  first  four  modes,  with  and  without 
piezoceramic  contribution  on  elements  1,  10,  11,  and  18.  Table  2  below  fists  the 
frequencies  associated  with  the  first  four  modes  generated  by  I-DEAS  FEA  for  the  arm 


Mode 

Frequency  (ELz) 
(without  piezos) 

Frequency  (Hz) 

(with  piezos) 

1 

0.136 

0.159 

2 

0.372 

0.469 

3 

1.651 

1.87 

4 

2.373 

3.06 

Table  2.  I-DEAS  FEA  Modal  Frequencies 


B.  MATLAB™  FEA  MODELING 

The  MATLAB  code  used  to  conduct  FEA  on  the  FSS  arm  appears  in  Appendix 
B.  No  control  laws  were  implemented  relating  sensor  output  to  actuator  input.  The 
code,  therefore,  performs  the  eigenvalue  problem  shown  in  Equation  (2).  This  code  was 
also  run  with  and  without  the  PZT  properties  included  in  order  to  compare  it  with  the 
I-DEAS  output.  Appendix  C  contains  printouts  of  the  modal  frequencies,  with  and 
without  piezos,  in  descending  order.  Table  3  below  fists  the  modal  frequency  predictions, 
for  the  first  four  modes,  generated  from  MATLAB.  The  I-DEAS  predictions  are  also 
re-listed  for  comparisons  sake. 
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Mode 

MATLAB 
(with  piezos) 

(Hz) 

I-DEAS 
(with  piezos) 

(Hz) 

MATLAB 

(without 

piezos) 

(Hz) 

I-DEAS 

(without 

piezos) 

(Hz) 

1 

0.162 

0.159 

0.136 

0.136 

2 

0.46 

0.469 

0.373 

0.372 

3 

1.787 

1.87 

1.646 

1.651 

4 

2.851 

3.06 

2.373 

2.373 

Table  3.  MATLAB  and  I-DEAS  FEA  Modal  Frequencies 


As  seen  from  the  table,  MATLAB  and  I-DEAS  modal  frequency  predictions  for 
both  cases  correlate  well.  For  the  case  with  piezo  contribution  to  the  arm  considered,  the 
I-DEAS  predictions  diverge  shghtly  from  MATLAB.  The  MATLAB  and  I-DEAS 
predictions  for  the  arm  without  piezo  contribution  considered  are  essentially  the  same. 
The  overall  effect  of  the  piezos  on  the  arm  stiffens  the  structure  incrementally,  as 
evidenced  by  the  shghtly  higher  frequencies  generated  by  both  MATLAB  and  I-DEAS. 
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IV.  MODAL  TESTING 


A.  INTRODUCTION 

Li  Chapters  n  and  HI,  the  differential  equations,  or  models,  of  the  system  are 
assumed  to  be  known,  and  the  theory  developed  consists  of  calculating  and  characterizing 
the  response  of  the  system  to  known  inputs.  This  is  called  the  "forward  problem".  In  this 
chapter  the  interest  hes  in  measuring  the  response  of  a  structure  and  in  some  way 
determining  the  equations  of  motion  from  the  test  data.  The  problem  of  determining  a 
system  of  equations  from  information  about  inputs  and  responses  belongs  to  a  class  of 
problems  called  "inverse  problems",  or  system  identification.  The  measurement  of  the 
vibrational  behavior  of  the  structure  is  used  to  verify  the  mathematical  model  of  the  test 
structure  outlined  in  Chapters  n  and  HI.  [Ref  3] 

Modal  testing  is  the  interpretation  of  test  data  collected  from  a  vibrating  structure. 
The  purpose  of  a  modal  test  is  to  construct  a  mathematical  model  of  the  vibrational 
properties  and  dynamic  behavior  of  the  tested  structure. 

Modal  testing  focuses  on  frequency,  shape  determination  and  damping.  On  a  plot 
of  frequency  response  as  a  fimction  of  frequency,  a  resonance  frequency  (modal 
frequency)  is  marked  by  a  peak  of  the  magnitude.  A  second  phenomena  of  resonance  is 
that  the  phase  of  the  response  shifts  by  180  degrees  as  the  frequency  sweeps  through 
resonance,  with  the  value  of  the  phase  at  resonance  being  90  degrees. 

B.  MEASUREMENT  HARDWARE 

Vibration  measurement  generally  requires  several  hardware  components.  These 
consist  of  a  source  of  excitation  for  providing  a  known  and  controlled  input  to  the 
structure,  a  transducer  or  accelerometer  which  converts  mechanical  motion  of  the 
structure  to  an  electrical  signal,  and  a  signal  conditioning  amplifier  to  match  the  signal 
level  of  the  accelerometer  to  that  of  the  data  acquisition  system  This  data  acquisition 
system  provides  the  interface  between  excitation  forces  and  accelerometer  signals  to  the 
testing  software  imbedded  within  I-DEAS. 
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Appendix  D  delineates  all  hardware  associated  with  modal  testing  of  the  FSS  arm. 
An  impact  hammer  with  a  built-in  force  transducer  was  used  to  excite  the  structure. 
Although  the  impact  hammer  is  simple  and  adds  no  mass  to  the  arm,  it  sometimes  may  be 
incapable  of  transforming  sufficient  energy  to  the  structure  to  develop  adequate  response 
signals  in  the  frequency  range  of  interest.  Consistent  impact  techniques  can  also  be 
difficult  to  develop.  The  key  frequency  range  of  interest  for  this  test  was  from  0  to  4  Hz. 
This  bandwidth  would  hopefully  contain  resonance  information  for  the  first  four  flexible 
modes. 

C.  DIGITAL  SIGNAL  PROCESSING 

Modal  Testing  analysis  is  performed  in  the  frequency  domain.  The  data  acquisition 
system  converts  the  analog  time  domain  signals  into  the  digital  frequency  domain  via  the 
Discrete  Fourier  transform.  This  data  is  then  passed  on  to  I-DEAS  in  order  to  perform 
the  required  computations.  Briefly,  by  comparing  the  transfer  functions  of  the  response  of 
the  structure  to  the  transfer  function  of  the  excitation  force,  one  can  determine  the  transfer 
function  of  the  structure  itself  and  ultimately  its  associated  resonances.  Mode  shapes  may 
also  be  extracted  from  this  data.  The  transfer  function  relationship  is  given  by; 

Gm  =  (23) 

where, 

Sfx(j&)=  response  of  the  structure  in  the  frequency  domain,  CO 
5^’to)=  excitation  force  in  the  frequency  domain,  0 
G{j(0)=  the  transfer  function  of  the  structure 

Reference  7  provides  an  excellent  source  for  information  on  the  theoretical  and  practical 
problems  surrounding  the  measurement,  processing,  and  analysis  of  modal  test  data. 


16 


D.  TEST  PROCEDURES 

A  modal  test  mesh  must  be  created  within  the  I-DEAS  testing  environment.  This 
mesh  is  depicted  in  Figure  6.  Unlike  the  mesh  created  for  the  finite  element  analysis,  no 
material  properties  need  be  entered.  The  nodal  coordinates  in  this  mesh  defines  the 
degrees  of  freedom  for  the  structure.  The  accelerometer  was  located  at  nodes  13,  10  and 
8  for  the  bulk  of  the  testing  due  to  physical  constraints  (proximity  to  the  piezos  and 
associated  wiring).  Node  16  was  used  as  the  primary  impact  location  for  defining  the  first 
four  modal  frequencies.  Impact  at  this  location  proved  to  consistently  and  reliably  excite 
the  first  four  modal  frequencies.  Impact  at  all  other  nodes  was  performed  in  order  to 
determine  mode  shapes.  [Ref  8] 
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Figure  6.  Modal  Test  Mesh  for 


Flexible  Spacecraft  Simulator  Arm 
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E.  TEST  RESULTS 


Figures  7,  8  and  9  are  I-DEAS  printouts  of  the  arm's  frequency  response  function 
measured  at  nodes  13,  10  and  8  respectively.  Impact  for  all  three  of  these  measurements 
took  place  at  node  16.  The  resonant  frequencies  are  clearly  represented  by  distinct  peaks. 
The  180  degree  phase  shift  (or  phase  reversal)  crossing  resonance  is  also  reasonably  well 
defiped.  Table  4  below  lists  the  averages  of  the  modal  frequencies  observed  during  these 
tests  at  nodes  13,  10  and  8.  The  table  also  re-hsts  the  MATLAB  FEA  predictions  (with 
piezoceramic  contribution)  and  a  percent  deviation. 


Mode 

Modal  Test 
Frequency  Averages 
(impact  at  node  16, 
response  at 
nodes  13,  10  and  8) 
(Hz) 

MATLAB  FEA 
Prediction 
(with  piezos) 

(Hz) 

% 

Deviation 

1 

0.124 

0.162 

23.4 

2 

0.436 

0.46 

5.2 

3 

1.374 

1.787 

22.8 

4 

2.751 

2.851 

3.5 

Table  4.  Modal  Test  Averages  and  %  Deviation  from  FEA 


Figure  10  depicts  an  analytical  cmvefit  that  was  performed  by  I-DEAS  on  the  test 
data  obtained  at  node  13.  This  analytical  curve  is  distinguished  from  the  raw  data  by  its 
smooth  character.  The  peaks  identified  by  this  curvefit  are  essentially  identical  to  those 
shown  on  Figure  7.  The  frequency  data  extracted  from  this  curvefit  was  input  to  the  mode 
shape  subroutine  within  the  1-DEAS  Testing  package. 

Figure  11  compares  the  third  mode  shape  generated  by  I-DEAS  FEA  to  that 
actually  observed  during  modal  testing.  The  bottom  illustration  of  Figure  11  is  the  third 
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Figure  9.  Modal  Test  Response  Function  (impact  at  node  16,  accelerometer  at  node  8) 


Figure  10.  Modal  Test  Response  Function  Analytical  Curvefit 
(impact  at  node  16,  accelerometer  at  node  13) 
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mode  shape  generated  during  the  testing,  the  top  represents  the  prediction  of  mode  3 
behavior  via  FEA. 

A  structure's  inherent  damping  characteristics  can  only  be  determined 

experimentally.  An  n  degree  of  freedom  system  has  n  damping  ratios, 
(dimensionless)  one  associated  with  each  modal  frequency,  CO,-,  (radians  per  second).  The 
free  vibration  characteristics  of  each  mode  decays  over  time,  t,  according  to  the  term: 

(23) 

[Ref  9].  Table  5  lists  the  damping  ratios  for  the  first  through  fourth  modes  obtained 
from  the  I-DEAS  testing  package.  The  exponential  time  decay  coefficient  is  also  listed. 


Mode,  i 

(0/  =  2n(freq) 
(radians/s) 

C/CO/ 

1 

0.779 

2 

2.739 

0.049 

3 

0.013 

8.633 

0.112 

4 

0.005 

17.285 

0.086 

Table  5.  Modal  Test  Damping  Ratios 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  major  reason  for  conducting  model  testing  is  for  comparison  and  verification 
of  analytical  models.  What  is  ultimately  desired  is  a  mathematical  model  of  the  structure 
under  consideration  for  the  purpose  of  predicting  how  the  structure  will  behave  under  a 
variety  of  different  loadings,  to  provide  the  plant  in  a  control  system  design,  such  as 
Positive  Position  Feedback  (PPF)  or  Strain  Rate  Feedback  (SRF),  and  to  aid  in  the  design 
process  in  general. 

Modal  testing  of  the  FSS  arm  provided  valuable  insight  into  this  field.  Testing 
results  yielded  resonant  frequencies  that  deviated  from  prediction  by  as  httle  as  3.5%  for 
mode  four,  and  as  much  as  23.4%  for  mode  one.  Overall,  modes  one  and  three  deviated 
from  prediction  on  the  average  of  23.1%,  and  modes  two  and  four  deviated  on  the 
average  of  4.4%.  The  difference  between  the  analytical  prediction  and  experimental 
determination  of  frequency  for  the  first  mode  is  significant.  Further  work  is  required  to 
improve  the  model  and/or  improve  the  modal  testing  technique. 

The  damping  ratios  obtained  during  modal  testmg  are  that  of  an  underdamped 
system,  (^,  <  1)  ,  [Refs  3  and  4].  The  expected  dominance  of  mode  one  is  supported 
by  the  lowest  exponential  time  decay  coefficient. 

As  mentioned  earher,  the  most  common  apphcation  of  modal  testing  is  the 
conparison  of  measured  vibration  modes  with  those  predicted  by  a  finite  element  or  other 
theoretical  models.  Resonance  frequencies  must  be  accurately  identified,  also  a  complete 
and  accurate  description  of  mode  shapes  is  desired,  especially  for  non-collocated  sensor 
and  actuator  systems.  Frequently,  the  analytical  model  will  not  predict  the  measured 
frequencies.  As  a  result  the  analytical  model  may  then  be  updated,  iteratively,  until  it 
produces  the  measured  natural  frequencies.  The  modified  model  is  then  considered  an 
inprovement  over  the  previous  model.  [Ref  3] 

One  possible  model  modification  would  be  to  consider  the  inertial  contributions  of 
the  mass  intensifiers  and  airpads.  These  were  treated  as  point  masses  with  zero  moment 
of  inertia.  This  would  be  a  good  starting  point  in  model  evolvement. 
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Normally,  the  FSS  arm  is  attached  to  the  center  body  of  the  FSS  which  houses  the 
reaction  wheel  and  thrusters  [Ref  2].  With  at  least  two  additional  "zero-Hz  operation" 
accelerometers,  such  as  that  used  for  this  experiment,  placed  on  the  arm  ( total  of  three 
on  the  arm),  and  with  an  additional  and  suitable  accelerometer  mounted  to  the  center  body 
for  excitation  source  measurement,  more  reliable  frequency  measurements  may  be 
expected.  These  results  may  then  be  used  to  "tweak"  any  implementation  of  PPF  or  SRF 
that  may  be  used  for  active  control,  in  more  of  a  real  time  fashion. 
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APPENDIX  A 
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without  piezos 


without  piezos 


without  piezos 


with  piezos 
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appendix  b 


%  FINITE  ELEMENT  ANALYSIS  OF  FSS  ARM 


clear 

m=l . 6375e-5;  %  Linear  Mass  Density  of  Beam  in  slugs/12/in 

mp=2.162e-5;  %  Linear  Mass  Density  of  Piezo  in  slugs/12/in 

%mp=0;  %  no  piezo  option 

h=2.5; 

hl=4.75; 

h8=4.25; 

hl2=4.8725; 

hl7=1.625; 

hl8=5.0; 

%  beam: 

I=2.0345e-5;  %  in"4 

E=1.0298e7;  %  lb/in''2 

%  piezo: 

Ip=4.56e-5;  %  in''4 

Ep=9.137e6/  %  lb/in''2 


%  Mass  Intensifiers  in  slugs/12: 

Lma=0 .00264 ; 

Lmb=0.0021; 

Lmc=0 .00472; 

Accel=0 .000416; 

n=18;  %  #  of  elements 


%  elemental  mass  and  stiffness 

mll=m*h*[156  22*h;22*h  4*h^2]/420; 
m22=m*h*[156  -22*h;-22*h  4*h''2]/420; 
ml2=m*h*[54  -13*h;13*h  -3*h"2]/420; 
m21=ml2 ' ; 


m22el=(m+mp) *hl* [156  -22*hl ; -22*hl  4*hl"2]/420; 
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mlle8=m*h8*  [156  22*h8;22*h8  4*h8''2] /420; 
m22e8=m*h8*  [156  -22*h8 ; -22*h8  4*h8''2] /420; 
ml2e8=m*h8* [54  -13*h8;13*h8  -3*h8-2]/420; 
m21e8=ml2e8'; 

mllell=(m+mp) *h* [156  22*h;22*h  4*h^2]/420; 
m22ell=(m+mp)  [156  -22*h;-22*h  4*h''2]/420; 
ml2ell=(m+mp) *h* [54  -13*h;13*h  -3*h^2]/420; 
ia21ell=ml2ell'  ; 

mllel2=m*hl2*  [156  22*hl2/ 22*hl2  4*hl2''2] /420; 
m22el2=m*hl2*  [156  -22*hl2 ; -22*hl2  4*hl2''2] /420; 
ml2el2=m*hl2* [54  -13*hl2; 13*hl2  -3*hl2^2] /420; 
m21el2=Inl2el2'  ; 

nillel7=m*hl7*  [156  22*hl7 ; 22*hl7  4*hl7^2] /420; 
Iii22el7=m*hl7*  [156  -22*hl7 ; -22*hl7  4*hl7''2] /420; 
ml2el7=m*hl7* [54  -13*hl7 ; 13*hl7  -3*hl7^2] /420; 
in21el7=ml2el7 '  ; 

mllel8=(iti+mp)  *hl8*  [156  22*hl8; 22*hl8  4*hl8"2]/420; 
m22el8=  (la+mp)  *hl8*  [  156  -22*hl8 ; -22*hl8  4*hl8^2] /420; 
ml2el8=(m+mp)*hl8*[54  -13*hl8 ;  13*hl8  -3*hl8''2 ] /420; 
m21el8=ml2el8'  ; 

mov=zeros (4)  ; 
mov ( 1 : 2 , 1 : 2 ) =m22 ; 
mov(2,2)=m22 (2, 2 ) +mllell (2,2) ; 
mov (2, 3)=ml2ell (2,1) ; 
mov (3, 2) =ml2ell (2,1)  ; 
mov (2,4) =ml2ell (2,2)  ; 
mov(4,2)=ml2ell (2,2) ; 
mov (3 : 4 , 3 ; 4 ) =m22ell ; 

kll=E*I*[12  6*h;6*h  4*h^2]/h"3; 
k22=E*I*[12  -6*h;-6*h  4*h-2]/h''3; 
kl2=E*I*[-12  6*h;-6*h  2*h''2 ] /h^3 ,• 
k21=kl2’  ; 

kpl=Ep*Ip*  [12  6*hl  -12  6*hl; 

6*hl  4*hl"2  -6*hl  2*hl''2; 

-12  -6*hl  12  -6*hl; 

6*hl  2*hl'^2  -6*hl  4*hl'^2  ] /hl'^3  ; 

kp2=Ep* Ip* [ 12  6*h  -12  6*h; 
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6*h  4*h^2  -6*h  2*h^2; 

-12  -6*h  12  -6*h; 

6*h  2*h^2  -6*h  4*h^2]/h^3; 


%  no  piezo  option 

%k;pl=zeros  (4)  ; 

%kp2=zeros (4) ; 

k22el=(E*I* [12  -6*hl;-6*hl  4*hl^2] /hl^3) +kpl (3 : 4, 3 : 4) ; 

klle8=E*I*  [12  6*h8;6*h8  4*h8^2] /h8''3; 
k22e8=E*I* [12  -6*h8;-6*h8  4*h8^2 ] /h8^3; 
kl2e8=E*I* [-12  6*h8;-6*h8  2*h8"2] /h8^3; 
k21e8=kl2e8  '  ; 

kllell=(E*I*  [12  6*h;6*h  4*h''2] /h"3) +kp2  (1 :  2, 1 ;  2)  ; 
k22ell=  (E*I*  [12  -6*h;-6*h  4*h^2] /h'^3) +kp2  (3  :  4, 3  :  4)  ; 
kl2ell=(E*I*  [-12  6*h;-6*h  2*h''2] /h^3) +kp2  (1 :  2, 3  :  4)  ; 
k21ell=kl2ell'  ; 

kllel2=E*I* [12  6*hl2;6*hl2  4*hl2^2] /hl2"3; 
k22el2=E*I*  [12  -6*hl2 ; -6*hl2  4*hl2^2] /hl2^3; 
kl2el2=E*I*  [-12  6*hl2;-6*hl2  2*hl2''2] /hl2''3; 
k21el2=kl2el2'  ; 

kllel7=E*I* [12  6*hl7;6*hl7  4*hl7^2] /hl7^3; 
k22el7=E*I* [12  -6*hl7 ; -6*hl7  4*hl7^2 ] /hl7^3; 
kl2el7=E*I* [-12  6*hl7;-6*hl7  2*hl7-2] /hl7-3; 
k21el7=kl2el7  '  ; 

kllel8=(E*I* [12  6*hl8;6*hl8  4*hl8^2] /hl8^3) +kpl ( 1 : 2, 1 : 2) ; 
k22el8=(E*I* [12  -6*hl8;-6*hl8  4*hl8^2] /hl8^3) +kpl (3: 4, 3: 4) ; 
kl2el8={E*I* [-12  6*hl8;-6*hl8  2*hl8-2] /hl8^3) +kpl (1 : 2, 3 : 4) ; 
k21el8=kl2el8' ; 

kov=zeros (4); 
kov(l:2,l:2)=k22; 
kov(2,2)=k22 (2,2)+kllell (2,2) ; 
kov(2,  3)=kl2ell  (2,1)  ; 
kov(3,2)=kl2ell (2,1) ; 
kov(2, 4)=kl2ell (2,2)  ; 
kov(4,2)=kl2ell (2,2) ; 
kov(3:4, 3;4)=k22ell; 
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%  global  mass  matrix 

gm ( 1 ; 2 , 1 : 2 ) =m2  2  e 1 +ml 1 ; 
gm ( 1 : 2, 3 : 4 ) =ml2; 

for  i=l:n-2 


gm(2*i+l:2*i+2,2*i-l:2*i)=m21; 

gm(2*i+l:2*i+2,2*i+l:2*i+2)=m22+mll; 

gm(2*i+l:2*i+2,2*i+3:2*i+4)=ml2; 


end 

gm(17 : 20, 17 : 20) =mov; 

gm(13: 14, 13:14) =m22+mlle8 ; 
gm(15:16, 13:14) =m21e8 ; 
gm(13:14,15:16) =ml2e8; 
gm(15:16, 15:16) =m22e8+mll ; 

gm  (19:20,  19:20) =m22ell+mllell ; 
gm(21:22, 19:20)=m21ell; 
gm(19:20, 21:22) =ml2ell ; 
gm  (21:22,21:22) =m2  2  e 1 1 +ml 1 e 1 2 ; 

gm (23:24,21:22) =m2 1 e 1 2 ; 
gm(21:22,23:24)=ml2el2; 
gm (23 : 24,23:24) =m22el2+mll ; 

gm ( 31 : 32 , 31 : 32 ) =m22+mllel7; 
gm(33:34,31:32)=m21el7; 
gm(31:32, 33:34 )=ml2el7 ; 
gm(33:34,33:34)=m22el7+mllel8; 

gm(35:36, 33:34) =m21el8 ; 
gm (33:34,35:36) =ml 2  e 1 8 ; 
gm (35:36,35:36) =m2  2  e 1 8 ; 

Lmm=zeros (size (gm) ) ; 

Lmm (3,3) =Lma ; 

Lmm (5,5) =Lma ; 

Lmm (7,7) =Lma ; 

Lmm (9, 9) =Lma; 
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FEA.M 


Lmm (11,11) =Lma ; 

Lram (13,13) =Lma / 

Lmm (23,23) =Lma ; 

Lmm (25,25) =Lma ; 

Lmm  (29,29)  =Lma ; 

Lmm (31,31) =Lma ; 

Lmm (27,27) =Accel ; 

Lmm (17, 17) =Lmb+0 .0154; 

Lmm (35, 35) =Lmc; 
gmm=gm+Lmm; 

%  global  stiffness  matrix 

gk(l:2, l:2)=k22el+kll; 
gk(l:2,3:4)=)cl2; 

for  i=l:n-2 

gk(2*i+l:2*i+2,2*i-l:2*i)=k21; 
gk(2*i+l:2*i+2,2*i+l:2*i+2)=k22+kll; 
gk(2*i+l:2*i+2, 2*i+3:2*i+4)=kl2; 

end 

gk(17:20, 17:20)=kov; 

gk(13: 14, 13 : 14 ) =k22+klle8 ; 
gk(15: 16, 13:14)=k21e8; 
gk(13:14,15:16)=kl2e8; 
gk(15:16, 15:16)=k22e8+kll; 

gk(19:20, 19:20)=k22ell+kllell; 
gk(21:22, 19:20)=k21ell; 
gk(19:20,21:22)=kl2ell; 
gk(21:22,21:22)=k22ell+kllel2; 

gk(23:24,21:22)=k21el2; 

gk(21:22,23:24)=kl2el2; 

gk(23:24,23:24)=k22el2+kll; 

gk(31:32,31:32)=k22+kllel7; 
gk(33:34,31:32)=k21el7; 
gk(31:32, 33:34)=kl2el7; 
gk(33:34,  33:34)=k22el7+kllel8; 

gk(35:36,33:34)=k21el8; 
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gk(33:34,35:36)=kl2el8; 

gk(35:36,35:36)=k22el8; 


%  mode  shapes  and  natural  freqs 

B=inv (gmm) *gk; 

[x,D]=eig  (B)  ; 
w=sqrt (D) ; 
f=w/2/pi; 

[kl, k2] =sort (diag (D) )  ; 
x=x ( : , k2 ) ; 

%  modal  masses 

y=x'*gmm*x; 

%  mode  shapes  normalized  wrt  mass 
for  i=l:2*n 

z ( : , i) =x ( : , i) /sqrt (y (i, i) )  ; 


end 

z; 
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»  f ea 
»  diag  (f) 


APPENDIX  c 


ans  =5 

l.Oe+003  * 

4.59118  474609514  Modal  Frequencies  (without  piezos),  descendins:  order 

3.99568135727380  r  o 

3.59480701359855 

2.99704503076688 

2.93678510562751 

2.28716854586410 

2.12001248476245 

2.00095892851113 

1.64377826366426 

1.57337814606062 

1.36667629704302 

1.24571931505966 

1.131393078.76285 

1.08173078771907 

0.83137661757685 

0.80808968513644 

0.67562468034393 

0.65960922618500 

0.33282677547923 

0.22223482048505 

0.206880508621L9 

0.18556683346412 

0.11445097218405 

0.09656670691312 

0.06764056549358 

0.04693032740926 

0.03658424100849 

0.02811923523049 

0.02120575875813 

0.01459733742682 

0.00694405426306 

0.00599862204817 

0.00237295384752 

0.00164563740851 

0.00037335596675 

0.00013563176832 
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»  diag(f) 
ans  = 


l.Oe+003  * 

4.48533942414282 

4.20598887360221 

3.93601883252608 

2.95655175921535 

2.91232463252813 

2.40665996231507 

2.12840833617064 

2.01547862614741 

1.67501528412810 

1.61012336453201 

1.39849731540630 

1.27385368688073 

1.14741402755860 

1.10615287290734 

0.85542324830144 

0.83293574716167 

0.72369166064969 

0.68994155774074 

0.34608426976888 

0.22304670115122 

0.22787352029346 

0.18559636091107 

0.11394978872364 

0.08854986624844 

0.06796197347123 

0.04772133056046 

0.03761481660823 

0.02913274059718 

0.02217253369177 

0.01562031009177 

0.00771814099641 

0.00652352892696 

0.00285149622307 

0.00178656363092 

0.00046041550163 

0.00016213258818 


Modal  Frequencies  (with  piezos),  descending  order 


» 
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APPENDIX  D.  MODAL  TESTING  HARDWARE 


A.  ACCELEROMETER 

1.  Model 
Sensotec  Model  #  3629-05 

2.  Available  Ranges 

+/-  lOg 

3.  Sensitivity 

10  mV/g 

4.  Usable  Frequency  (Bandwidth) 

0-250  Hz 

5.  Excitation 

+/-5V 

B.  SIGNAL  CONDITIONER 

1.  Model 

Instruments  Division  Measurement  Group,  Inc.  Model  2100 

2.  Type 

Strain  Gage  Conditioner  and  Amplifier  System  with  Integral  Power  Supply 

3.  Gain 
1-2100 

4.  Operation 

Full  Balanced  Operation  for  Noise  Rejection 

C.  IMPULSE  FORCE  HAMMER 

1.  Model 

PCB  Piezotronics  Model  #  086B01 

2.  Range 
0-100  Lb. 
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3.  Hammer  Sensitivity 
52.2  mV/g 

D.  HAMMER  POWER  SUPPLY 

PCB  Piezotronics  Model  #  483B03  Line  Power  Supply 

E.  DATA  ACQUISITION  HARDWARE 

Hewlett-Packard  Model  HP  3  5  65  S  Measurement/Data  Acquisition  System 
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